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Abstract: We discuss the problem of theoretical uncertainties in the combination 
of observables related to the CKM matrix elements and propose a statistically sen- 
sible method for combining them. The overall fit is performed on present data, 



and constraints on the matrix elements are presented as well as on /^^yS^^. We 
then explore the implications of recent measurements and developments: J/ip CP 
asymmetry, e'/e and B Kir branching fractions. Finally, we extract from the 
overall fit the Standard Model expectations for the rare kaon decays K — > vrz/z/. 
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1. Introduction 

The Cabbibo-Kobayashi-Maskawa (CKM) matrix is extensively studied nowadays. 
With the birth of the new B factories and the upgrade of the Tevatron experiments, 
high precision measurements in the B meson sector are expected and the question 
of testing the CKM ansatz is pushed towards more and more stringent hmits. By 
"testing" we mean two aspects: 

1. Given that the Standard Model (SM) is right, what is the best knowledge we 
have on the CKM free parameters? 

2. Are all the measurements involving CKM matrix elements compatible within 
their errors? 

As we shall see (section 2.2), a key point in performing this test is a proper treatment 
of the theoretical estimates that enter the description of the observables. 

First attempts to combine several observables were performed by simply drawing 
in the Unitarity Triangle (UT) plane individual 95% CL regions for p and obtained 
by varying coherently the experimental and theoretical errors of each observable. 
These regions had the advantage of being geometrically simple^ . The intersection of 
these regions was taken as a 95% CL for (p, rf) . While statistically wrong (it neglects 
the correlations induced by the combination), this method gives surprisingly good 
results. It can certainly be used to get an insight into the effect on CKM parameters 
of a given set of observables. 

More sophisticated fits have been proposed [1, 2] but, as in the previous case, 
assuming some flat (or Gaussian) distribution for theoretical parameters. 

We propose here a way of decoupling the experimental measurements from the 
theoretical estimates (section 2.3). It will be illustrated in the (^,77) plane and the 
(sin2Q;, sin 2/5) plane (section 3). The overall fit leads also to constraints on the 
theoretical parameters fBd^Bs^ and (section 3.4). Recent experimental and 
theoretical developments are then investigated (section 4). Finally, the impact of B 
factories is discussed and the SM rare K branching ratios are extracted (section 5). 

2. Outline of the method 
2.1 Least Squares method 

The method exposed here is discussed in more detail in [3]. It can naturally accom- 
modate any new measurements (as in the case of the Fleischer-Mannel bound [4]). 

The CKM matrix contains 4 independent parameters, which can be taken as 
three angles and a phase, or, exhibiting a hierarchy, as the (improved) Wolfenstein 

^Note that this is no longer the case when considering for instance sin 2a. 
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parameters: A, A, p, fj. The Cabbibo angle is already very well known [5] and A = 
sin 6c = 0.2205 will be considered as fixed in the following. 

Our taslc is therefore to determine 3 independent parameters, given a set of 
measurements Fj =< Fj > icTj and their theoretical description Y{A,p,fi). 

Considering the two aspects discussed in the introduction, we are naturally lead 
to use the least squares estimate method. One builds: 



<Yi>-Y{A,p,fj) 



(2.1) 



Minimising the function: 



1. the values of {A,p,f]) taken at the minimum of the function provide the least 
square estimates. Hyper-regions (as the {p,fj) projection, i.e. the UT) at a 
given confidence level can be constructed. 

2. To test the compatibility between all measurements (and their theoretical de- 
scription), one studies the value taken by the function at its minimum, Xmin- 
When the errors are Gaussian, one can further use the probability distribu- 
tion to quantify it. 

2.2 The problem with theoretical estimates 

Life would be simple if factors that include some level of model dependence did not 
enter into the calculation. In the following, these will the bag factors for K and 

B mesons {Bk,Bb)- the B decay constant (/s, which will be combined with the 
previous bag factor to give fBi\jBB^) and a large part of the error related to \Vub\- 

While the error on a measured quantity has a clear statistical meaning^, this 
is not the case for the errors quoted by theorists for their estimates. Here enters a 
level of subjectivity reflecting a "degree of belief" that the true value lies inside some 
range [6]. This however may vary from one person to another and is impossible to 
quantify in terms of a probability distribution function (p.d.f.). 

Furthermore, not knowing what the p.d.f. of a given theoretical estimate is makes 
it impossible to combine with the other estimates on the market. 

Finally, a crucial point is that, not knowing where the true value of a parameter 
lies is not equivalent to taking a flat probability distribution within some bounds. This 
is not a matter of philosophical discussion; any p.d.f. extracted from a combination 
using a fiat p. d.f for theoretical estimates is senseless. 

2.3 One solution 

One method to overcome these problems is the following: 



^To some level one may argue that systematic errors have also unknown p.d.f. - whenever possible 
we will therefore include this part of the error in the "model-dependent" term. 
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• Determine a "reasonable range" for each theoretical parameter, given the spread 
of published results. Here reasonable might mean "conservative" . 

• Bin the whole range, and scan all the values of the theoretical parameters. For 
each set of values (a "model"), build and minimise the ^ (equation 2.1) using 
statistical errors only. Since these are Gaussian, one can test the compatibility 
between the measurements /or this set of theoretical values. The estimate is 
rejected if it does not pass a Prob{x^) test. If it succeeds, draw a 95% CL 
contour in the UT plane. 

• Now, one does not know which contour is right, but presumably one of them 
will be correct (assuming the scanned range is really reasonable). Therefore 
our maximum knowledge is that the set of all these contours is an overall 95% 
CL. Visually, it corresponds to the envelope of all the individual contours. 

The only output of such a method is to provide an overall 95% CL region in the (p, f]) 
plane (or in (sin2Q;, sin 2/?)). No particular p.d.f. can be inferred, and estimates of 
the mean and the "cr" have no particular meaning; the p.d.f. within the envelope is 
simply unknown. 



3. CKM 1999 
3.1 Observables 

Our present knowledge of the CKM elements hes in the following observables. 

3.1.1 \Vcb\ 

\Vcb\ gives a direct access to the A parameter of the Wolfenstein parametrization. 
Much work on the experimental/theoretical side allows us to quote [8]: 

I K&l = 0.040 ±0.002 (3.1) 

Note that part of the error quoted here is indeed model-dependent, but it was checked 
that no "visual" difference can be observed when treating it as a model-dependent 
range. We can safely consider it as statistical in the following. 

3.1.2 \Vub\ 

• there has been a recent update of the CLEO B — > piu analysis [7] giving: 

iKbl = [3.25 ± O.U{stat)1:oll{syst) 



±0.55{model)] x 10 



-3 
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• also LEP reports a combined value [8] : 

iKfcl = [4.05t'ofdstat + rfet)l°f (6 ^ c) 
-oirib -^u)± 0.16{HQE)] x lO'^ 

where the last three terms corresponds to systematics related to the modelling 
of b ^ c, b ^ u and to the Heavy Quark expansion, used to relate the measured 
BR{b XJv) to \Vub\. 



The LEP measurements have a statistical error about twice as large as CLEO, 
and both experiments have the same order of magnitude for detector systematics. 
But clearly the dominant part of the error is model-dependent. Since it would be 
questionable to go beyond 15% (relative) error for this term [6] we will use as a 
reasonable guess estimate^: 

\Vub\ =< \Vub\ > ±0.13 X 10-^{stat) 

where we will vary the mean within the range: 

< iKbl >e [2.9,3.9] X 10"^ (3.2) 

Note that the final results do not depend crucially on the details of the value used 
here. 

3.1.3 I^k] 

We use for this measurement the PDG value [5] : 

le^l = (2.285 ± 0.018) x 10"^ (3.3) 

The theoretical description of this observable in terms of CKM matrix elements can 
be found in the literature (see for instance [9]) and requires a value for the non- 
perturbative QCD bag factor, Bk- Following Buras [9], we will use: 

Bk e [0.65, 0.95] (3.4) 

3.1.4 AruB, 

The LEP oscillation working group has produced an optimal combined value for the 
B mixing frequency [10]: 

AruBa = 0.473 ± 0.016 ps~^ (3.5) 

^We do not use in this part the CLEO inclusive analysis since it gives comparable results and is 
older. 
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Again, theoretical computations require not only a value for the corresponding bag 
factor Bb, but also in this case the B decay constant fs, so that finally the relevant 



model-dependent theoretical parameter is fs^sj Bb^- 

Here again, opinions may vary on a conservative range. Following the work 
performed in the BaBar Physics book [3], we will use: 



/fi.V^e [160,240]MeV 



(3.6) 



Given the large variations observed with the first unquenched estimates [11], this 
range may even be over-optimistic. 

3.1.5 Am^^ 

The study of the mixing frequency in the strange B meson sector has not lead to a 
firm measurement, but much information can be inferred from the combination of 
the amplitudes performed by the LEP Oscillation working group [10]. 

As already detailed in [3] and [4], we use the available information optimally by 
building a of the form: 



A{/\mBXAp,f])) 



(3.7) 



where A, are extracted from the amplitude curve, and Ams^iA, p,fj) is the the- 
oretical computation. 

It is frequently argued that theoretical errors cancel when taking the ratio 
AtubJ Atub^- This is only true to the level of precision determined by the new 
theoretical parameter 



(3.8) 



For this parameter, we will scan the range^ ([9] and references therein): 

€[1.12,1.48] 



(3.9) 



3.2 Unitarity Triangle 

We build the x^i^^P^v) defined in section 2.1 with all the observables described 
in the previous part, only using statistical errors (table 1). Each model-dependent 
parameter is scanned independently within the range of table 2. For each set of these 
values, the is minimised using the package MINUIT [12] and the estimate is kept 

^Note that since we basically have a lower bound on Attib^, the only relevant value for the 
CKM combination is the upper value £,gMAX ■ This factor may vary depending on the authors and, 
starting from similar "guesstimates" , the square factor enhances the differences. 
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Measurement 


Mean value 


Error 


\Vcb\ 


.040 


0.002 


IKbKx 103) 


< \Vuh\ > 


0.13 


\eK\{x 10^) 


2.285 


0.018 


Ams,(ps-i) 


0.473 


0.016 


AmB,c/.(3.1.5) 


A 





Table 1: The set of measured values used in the global fit. Whenever possible the error 
is statistical only. 



Parameter 


Min. 


Max. 


< IKfel > (xio3) 


2.9 


3.9 


Bk 


0.65 


0.95 




160 


240 


es 


1.12 


1.48 



Table 2: Range scanned in the global fit of the model-dependent theoretical parameters. 

if it satisfies a probability cut, -P(Xmm) ^ 0-05. We then draw the associated 
95% CL contour in the (p, f]) plane, for this model. Each model surviving the cut is 
superimposed. The envelope of all the contours is the overall 95% CL region for the 
CKM parameters {p^fj). 

From figure 1 one can extract (roughly) the projections: 

< p < 0.3 (3.10) 
0.2 < ^ < 0.45 (3.11) 

Prom such a global fit, an estimate of the third CKM parameter A is possible. 
However, given the model dependency induced on that parameter by the various 
observables, it is certainly wiser to extract A directly from the \Vcb\ measurement 
alone: 

A = = 0.82 ± 0.04 (3.12) 
A 

3.3 (sin2Q;, sin 2/3) 

The same can be built in another basis, namely (^4, sin 2q;, sin 2/3). The same 
procedure is apphed^ and figure 2 shows the 95% CL region in the (sin 2q;, sin 2/3) 
plane. 

^For completeness, there appears a four-fold ambiguity which is solved by taking the minimum 
value of the under the four hypotheses for each point [3]. 
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Prom that figure we get the projections (95% CL): 

.50 < sin 2/3 < .85 (3.13) 
-.95 < sin2Q; < .50 (3.14) 

The 95% CL regions we obtain are larger than those reported in [2], especially for the 
sin 2q parameter. This comes from a (somewhat) different choice of the parameters 
(mainly CsMAx) but especially from a different treatment of the theoretical errors; in 
[2] the 95% CL region for sin2Q; is extracted from the "p.d.f." inferred from the fit. 
As detailed in section 2.2, we disagree with that approach. 

3.4 Constraints on /^^ ^ Bb^ 

So far, we have just explored the first aspect of testing (getting the best knowledge 
on the CKM parameters assuming the SM is right) and we turn now to the second; 
are the results consistent? 

Here, recall the procedure; model dependent terms are scanned within a range 
and for each set of them, the is computed. The value of the at its minimum 
indicates the consistency of all measurements for that set of theoretical parameters. 
One can therefore reject sets of values which are inconsistent with all the measure- 
ments (if they were all rejected, we would conclude there is a consistency problem, 
implying new physics). 

Figure 3 shows a projection of all the scanned points that survived the Prob{xfnin) 
cut of 5% during the combination in the {fB^\jBB^-, Bk) plane. (Note these param- 
eters are somewhat related by lattice computations.) 

Figure 3 indicates that the Standard Model combination imphes: 

• /b,/b^> 195MeV 

• low values of with large values of fsd^jBEa are disfavoured 

4. Recent developments 

4.1 a^Ks (CDF) 

CDF has reported a first measurement of the J/il) CP asymmetry, which leads to 
the (model- independent) measurement [13]: 

sin2/3 = 0.79l°:^l (4.1) 

With respect to what is already known on sin 2/3 from the combination above (fig- 
ure 2), it clearly does not constrain the CKM matrix elements any further. 

This is however the first measurement directly related to the phases of the CKM 
matrix only and the fact that sin 2/3 > (95% CL) is strong support for the vahdity 
of the CKM description. 
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4.2 e'je (KTeV, NA48) 

A large value for this observable has been measured by the KTeV and NA48 collab- 
orations. Their combination gives [14] 

= (21.3±2.8) X 10-^ (4.2) 

The CKM description of this observable involves only r\: 

e'/e = r]A^X^F,, (4.3) 

but the function F^/ includes many theoretical factors. A crude description is [15]: 

'no MeV' 



= 13 



^ ms{mc) J ^4UiViev 

2.5 



340 MeV 



x[BW^)(l-n„^)-0.4Sf^'(^^|5gij 1 (4.4) 

The main uncertainties in this formula are related to the QCD penguins (Sg^^^^), 
electroweak penguins (Sg^^^^) and the strange quark mass (ms(mc)). 

The difficulty in predicting any value for this parameter comes from the fact 

('1/2) ('3/2) 

that two badly known parameters {B\ and ) are subtracted and that the 

difference could be as low as 0. On the other hand, this difference cannot be too large. 
And an upper bound on F^i directly translates into a lower bound on r] (equation 4.3). 
For instance using [15]: B^^'"^^ = 1.0±0.3, = 0.8±0.2, m,(mc) = 130±25MeV, 

a5^^^ = 340 ± 50 MeV, = 0.25 ± 0.08, rut = 165 ± 5 GeV and coherently varying 
the errors, one obtains: 



F^, < 13.4 (4.5) 



which translates into: 



and using the experimental input (sections 3.1 and 4.2), allowing a 2a variation, one 
obtains: 

77 > 0.32 (4.7) 

Given figure 1 this would be a strong constraint on the UT. 

Unfortunately, the formula (equation 4.4) is not accurate enough. But the mes- 
sage is; since the measurements provide a large value of e'/e, a firm upper bound on 
the theoretical parameters can be enough to constrain significantly the 77 parameter 
of the CKM matrix. 
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4.3 Bounds from B —>■ Kn 

4.3.1 The "Fleischer- Mannel bound" 

Fleischer and Mannel have proposed [16] a constraint on the angle 7 of the UT, solely 
from the measurements of CP- averaged branching ratio B — > Ktt: 

sin^ 7 < i? (4.8) 

where 

n^m^iiii^ (49) 

While there has been a lot of discussion about possible theoretical uncertainties that 
would weaken this bound [17] , it suffers mainly from the present CLEO measurement 
which gives a value of R consistent with one [19]: 

R = l-lliasi (4.10) 
It will certainly worth revisiting it when more accurate measurements become avail- 
able. 

4.3.2 The "Neubert-Rosner bound" 

Following that idea, Neubert and Rosner have proposed a bound on 7 from charged 
B decays only, in which theoretical uncertainties are much more under control [18]. 
It relies on two absolute branching ratios, which are used to define: 

r{B^^n^ 
^* - 2(S± ^ Tro^i) ^^-^^^ 

and 

" T{B^ ^ tt^tt") " ' 
r(5± ^ 7r±iro) 

where in this second formula, 9c is the Cabbibo angle and Rsuis) is a precisely known 
correction [20]. Using these observables the bound reads [21]: 

\fR^-l 

\Xr\ = \^ 1 < \dEw - C0S7I (4.13) 

£3/2 

where 5ew is calculable in terms of Standard Model parameters: 

Sew = (0.64 ± 0.09) x (4.14) 

I \Vcb\ I 

The bound in equation 4.13 is discriminant provided that Xr is (statistically) below 
one. While the first CLEO results were promising, the latest updates [19] do not 
confirm a value statistically different from one: 

Xr = 0.72 ± 0.98(exp) ± 0.03(th) (4.15) 

As in the Fleischer-Mannel case, one waits eagerly for more precise measurements 
from CLEO-III, Belle and BaBar. 



€3/2 = y2i?s{/(3) tan 6'c 



(4.12) 
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5. Future measurements 

5.1 B factories 

With the start of B factories, one can expect some new measurements of CP asym- 
metries which are related to sin 2/3 and sin2Q;. The stakes are however different. 
Given figure 2: 

• sin 2(3 is already well constrained. The goal of measuring it is to test the SM, 
since in a variety of models new physics may appear only in the CKM phases 
[22]. It is not expected that measuring this angle will constrain fBa\J much 
more than presently. 

• sin2Q; is largely unknown and the goal of B factories is to measure it. The 
extraction of that angle from the measured CP asymmetries is difficult (im- 
possible?) and will certainly require several years of running [3]. The most 
promising channel is presently i? — Svr in which all the amplitudes can be 
extracted from a global fit to the Dalitz plot. 

5.2 K TTvP 

Both the charged mode and the neutral one are "theoretically clean" and measuring 
their rate would significantly constrain the CKM matrix [9]. For the time being, 
we can extract the expected branching ratio from the global fit by scanning all the 
points in the contours of figure 1 and keeping the minimum and maximum value of 
the corresponding computed branching ratio. One obtains (for 95% CL): 

BR{Kl n^vv) G [1 - 5] X 10"^^ (5.1) 
BR{K+ Ti+uv) e [4 - 10] X 10"^^ (5.2) 

6. Conclusions 

We want to draw the attention of the reader to the difficulties that arise when includ- 
ing some theoretical estimates into an overall combination of observables relevant to 
the CKM determination. This was also emphasised by Stone [6] and will (and already 
does) limit our understanding of the CKM parameters. Given that a "theoretical" 
error has an unclear statistical meaning, we conclude that extracting any p.d.f. from 
a combination including these parameters is just meaningless, and that no "central 
values" and "errors" should be ever quoted. 

Nevertheless, we have proposed a (conservative) method to obtain some 95% CL 
regions for all CKM parameters, by separating the statistical errors due to measure- 
ments from the systematic and model-dependent ones. Prom such a combination, we 
obtain the 95% CL bounds: 
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• < p < 0.3, 0.2 < ?7 < 0.45 



• .50 < sin 2/3 < .85, - .95 < sm2Q; < .50 

• /b,/b^> 195MeV 

Among new developments, the large value measured for e' je could constrain 77 
(by a lower bound) if theoretical uncertainties were more under control (an upper 
bound would be sufficient). 

B K-K absolute branching ratios could constrain significantly the angle 7 of 
the UT but must be measured more precisely. 

Finally, from the overall combination, one can extract the expected branching 
ratios: 

BR{Kl ttVp) e [1 - 5] X 10"^^ 
BR{K+ n+uP) e [4 - 10] x 10"^^ 
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Figure 1: The allowed 95% CL region for p-ry. Each contour corresponds to one theoretical 
model, and the envelope of all of them is the overall 95% CL combination. Also shown (for 
historical reasons) the usual individual constraints. 
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Figure 2: The 95% CL allowed region for sin2a-sin2/3. Each contour corresponds to one 
theoretical model, and the envelope of all of them is the overall 95% CL combination. 
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Figure 3: The projection on tiie /b^ ^JBb^, Bk plane of the theoretical models that survive 
the consistency cut of 5% (points). The grey area depicts the whole range scanned. 



16 



